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Abstract. A fractafold, a space that is locally modeled on a specified fractal, is the fractal 
equivalent of a manifold. For compact fractafolds based on the Sierpihski gasket, it was 
shown by the first author how to compute the discrete spectrum of the Laplacian in terms 
of the spectrum of a finite graph Laplacian. A similar problem was solved by the second 
author for the case of infinite blowups of a Sierpihski gasket, where spectrum is pure point 
of infinite multiplicity. Both works used the method of spectral decimations to obtain 
explicit description of the eigenvalues and eigenfunctions. In this paper we combine the ideas 
from these earlier works to obtain a description of the spectral resolution of the Laplacian 
for noncompact fractafolds. Our main abstract results enable us to obtain a completely 
explicit description of the spectral resolution of the fractafold Laplacian. For some specific 
examples we turn the spectral resolution into a "Plancherel formula" . We also present such 
a formula for the graph Laplacian on the 3-regular tree, which appears to be a new result 
of independent interest. In the end we discuss periodic fractafolds and fractal fields. 
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1. Introduction 

Analysis on fractals has been developed based on the construction of Laplacians on certain 
basic fractals, such as the Sierpihski gasket, the Vicsek set, the Sierpihski carpet, etc., which 
may be regarded as generalizations of the unit interval, in that they are both compact and 
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have nonempty boundary. As is well-known in classical analysis, it is often more interesting 
and sometimes simpler to deal with spaces like the circle and the line, which have no bound- 
ary, and need not be compact. The theory of analysis on manifolds is the natural context 
for such investigations. The notion of fractafold, introduced in [57], is simply the fractal 
equivalent: a space that is locally modeled on a specified fractal. For compact fractafolds 
based on the Sierpihski gasket, it was shown in [37J how to compute the spectrum of the 
Laplacian in terms of the spectrum of a Laplacian on a graph T that describes how copies of 
SG are glued together to make the fractafold. On the other hand, in [JT] a similar problem 
was solved for the case of infinite blowups of SG. These are noncompact fractafolds where 
the graph T mirrors the self-similar structure of SG. Not surprisingly, the spectrum in the 
compact case is discrete, and the eigenvalues and eigenf unctions are described by the method 
of spectral decimation introduced in [TTJ. The surprise is that for the infinite blowups the 
spectrum is pure point, meaning that there is a basis of L 2 eigenf unctions (in fact compactly 
supported), but each eigenspace is infinite dimensional and the closure of the set of eigenval- 
ues is a Cantor set. Again the method of spectral decimations allows an explicit description 
of the eigenvalues and eigenfunctions. 

In this paper we combine the ideas from these earlier works [37| W\\ to obtain a description 
of the spectral resolution of the Laplacian for noncompact fractafolds with infinite cell graphs 
T. The graph Y is assumed to be 3-regular, so the fractafold has no boundary. The edge 
graph Tq is then 4-regular, and the fractafold is obtained as a limit of graphs obtained 
inductively from T Q by filling in detail (that is, each graph triangle is eventually replaced 
with a copy of the Sierpihski gasket). Our first main abstract result is Theorem 12. 3[ which 
describes how to obtain the spectral resolution of the Laplacian on the fractafold from the 
spectral resolution of the graph Laplacian on r . This is a version of spectral decimation, 
and uses an idea from [27J to control the L 2 norms of functions under spectral decimation. 
The second main abstract result is Theorem I3.1[ which shows how to obtain the spectral 
resolution of the graph Laplacian on To from the spectral resolution of the graph Laplacian 
on r using ideas from [M| HQ] . We note that the spectral resolution on To may or may not 
contain the discrete eigenvalues equal to 6, and the explicit determination of the 6-eigenspace 
and its eigenprojector must be determined in a case-by-case manner. Combining the two 
theorems enables us to obtain a completely explicit description of the spectral resolution of 
the fractafold Laplacian to the extent that we are able to solve the following problems: 

(a) Find the explicit spectral resolution of the graph Laplacian on T; 

(b) Find an explicit description of the 6-eigenspace and its eigenprojector for the graph 
Laplacian on To- 

The bulk of this paper is devoted to solving these two problems for some specific examples. 
However, we would like to highlight another problem that arises if we wish to turn a spectral 
resolution into a "Plancherel formula" . Typically we will write our spectral resolutions as 



where P(A, x, y) is an explicit kernel realizing the projection onto the A-eigenspace, i.e. 
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and dm(X) is a scalar spectral measure. (Here neither P(X,x,y) nor dm(X) are uniquely 
determined, since we can clearly multiply them by reciprocal functions of A while preserving 
(OJ) and ([L2l).) If we write 

(1.3) Pxf(x) = f P(X,x,y)f(y)dv(y) 

then (II .ip resolves / into its components P\f in the A-eigenspaces. A Plancherel formula 
would express the squared L 2 -norm ||/||| in terms of an integral of contributions from the 
components P\f. In the case of pure point spectrum this is straightforward, for then the 
A-integral is a discrete possibly infinite sum, and we just have to take the L 2 -norm of each 
Pa/, so 

(1-4) ||/||1= £ W P M 

Ag<x(-A) 

where Pa is the eigenprojection. The spectral measure m is the counting measure in this 
case. 

In the case of a continuous spectrum this is decidedly not correct, and there does not 
appear to be a generic method to obtain the correct analog. So we pose this as a third 
problem: 

(c) describe explicitly a Hilbert space of A-eigenfunctions with norm || ||a such that 
||Pa/||a is finite for m-a.e. A and 

(1.5) 11/112=/ \\Pxf\\ldm(X). 

This problem is interesting essentially only when the eigenspace is infinite dimensional. The 
resolution of this problem in some classical settings is discussed in [35J and [Hj. Here we 
present a solution to this problem for the graph Laplacian on the 3-regular tree. This result 
appears to be new, and is of independent interest. 

The first specific examples we consider is the tree fractafold, discussed in Section HJ where 
T is the 3-regular tree. In this case the solution to a) is well-known [U We solve (b) 
by showing that the 6-eigenspace on T is infinite dimensional and we give an explicit tight 
frame for this space. We solve (c) in terms of a mean value on the tree that is in fact different 
from the obvious mean value. The fractafold spectrum in this example is a union of point 
spectrum and absolutely continuous spectrum. 

In Section[5]we discuss periodic fractafolds, concentrating on a honeycomb fractafold where 
T is a hexagonal lattice. In this case the solution to a) is also well-known. Our solution to 
b) gives a basis for the infinite dimensional 6-eigenspace of compactly supported functions. 
Finally in Section [HI we discuss an example of a finitely ramified periodic Sierpihski fractal 
field (see [12]) that is not a fractafold, but can be treated using our methods. 

Essentially all the results of this paper can be extended to fractafolds based on the n- 
dimensional Sierpihski gasket, using similar methods. It seems likely that similar results 
could be obtained for any p.c.f. fractal for which the method of spectral decimation applies 
(see [H El [101 [TH H3l H6l UHl [SOI |26l EH [39l S2l S3! and references therein] ) . 

Acknowledgments. We are grateful to Peter Kuchment and Daniel Lenz for very help- 
ful discussions, and to Eugene B. Dynkin for asking questions about the periodic fractal 
structures. We thank Matthew Begue for help in the manuscript preparation. 
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2. Set-up results for infinite Sierpinski fractafolds 

2.1. Laplacian on the Sierpinski gasket. We denote by Asg the standard Laplacian 
on SG, and by hsg the standard normalized Hausdorff probability measure on SG (see 
[T71EH1ES] for details). The Laplacian A S g is self-adjoint on L 2 (SG, Hsg) with appropriate 




Figure 2.1. Sierpinski gasket. 

boundary conditions (usually Dirichlet or Neumann). The Laplacian A$g can be defined 
either probabilistically or analytically, using Kigami's resistance (or energy) form and the 
relation 



£(/,/) = -§ / f^SGfdpiSG 
1 JSG 



'SG 

for functions in the corresponding domain of the Laplacian. The energy is defined by 



£(/,/) 



lim (Tf E (m-f(y)f- 



In these formulas V n is a finite set of (3 n+1 + 3)/2 points in SG that are at the Euclidean dis- 
tance 2~ n from the neighboring points, and ~ denotes the recursively defined graph structure 
on V n . Note the normalization factor | that is inserted here for the convenience of compu- 
tation (see [39] the explanations). 

2.2. Spectral decimation and the eigenfunction extension map. Both Dirichlet and 
Neumann spectra of Asg are well known (see [Hi [39j SI]). To compute the spectrum of 
Asg one employs the so called spectral decimation method using inverse iterations of the 
polynomial 

R(z) = z{5-z). 

By convention the eigenvalue equation is written — A S gu = Xu because — Asg is a non- 
negative operator. Each positive eigenvalue can be written as 



(2.1) A = lim 5 m X m = 5 m ° lim 5 fc A 



k— >oo 



fc+mo 



for a sequence {A m }^ = such that A m = R(X m+ i) and X mo G {2, 5, 6}, which can be written 



I m=mo 

as 

jofc 



R ok (X k+mo ) G {2,5,6}, 

where the powers R ok of R are composition powers. If we denote %Kk{z) — R ok {5~ h z) then 



(2.2) R ok (X k+mo ) = D\ k (5 k X k+mo ) = ft k -5- mo 5 mo 5 fc A 



2 

- — iih. i r 1 1 'i ) r r< \ 



,3 

Thus an important role is played by the function 

(2.3) Vt(z) = lim R ok (5- k z). 

This is an analytic function, which is a classical object in complex dynamics, and a recent 
detailed study and background can be found in [6j [7j. In the context of the Laplacian on the 
Sierpinski gasket this function first appeared in [281 Lemma 2.1] and [SI Remark 2.5] (see 
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also [151 129] for related results). In particular, this function can be defined as the solution 
of the classical functional equation 

(2.4) R(^(z)) = m(5z). 

Note that, in a neighborhood of zero, the inverse of the function £R can be defined by 

(2.5) R(w) = lim 5 k R~ k (w), 

k— >oo 

and satisfies the functional equations 

(2.6) 5%(w) = R{R{w)), 

in a neighborhood of zero. 

One can see from ( 12. 2 p that each nonzero eigenvalue A satisfies 



A G 5 m °9T 1 {2, 5, 6} C |J 5 m 9r 1 {2,5,6}. 



m=0 



Some of the points in the union of these sets are so-called "forbidden eigenvalues" , and the 
rest are so-called 2-series, 5-series and 6-series eigenvalues (see [39] ) . A detailed analysis 
shows that the spectrum of the Dirichlet Laplacian is 

(oo 
mr 1 {2, 5} U 5JK -1 {5} |J S^fH' 1 ^^} 
m =2 

and the spectrum of the Neumann Laplacian is 

Sjv = {0}U5 (*H- 1 {3}U (J 5 m °9T 1 {3,5} J . 

\ mo=l / 

The multiplicities, which grow exponentially fast with k, were computed explicitly in [llj . 
and also can be found in [TJ EH E]. Note that, because of the functional equations (12.41) 
and (12.61) . and because R(2) = R(3) = 6, we have 

5 (m- 1 {2} U = ^ _1 {6}. 

If we define 

£ e ** = 5 |^- 1 {2}U Q 5 m *H~ 1 {5} J cST^e}. 



m=0 



then we have the following proposition. 



Proposition 2.1. For any v G dSG and any complex number A ^ Y> ext there is a unique 
continuous function 4>v,\(') '■ SG — > M, such that ip Vy \{v) = 1, ipv,\ vanishes at the other two 
boundary points, and the pointwise eigenfunction equation —Aip Vt \(x) = \ip Vt \(x) holds at 
every point x G SG\dSG. 

Naturally, ip Vj \ is called the eigenfunction extension map, which is explained in [23 Section 
3.2], and the proposition is essentially the same as [3"9| Theorem 3.2.2]. 

Example 2.2. Spectral decimation for the unit interval [0,1]. In order to illustrate these 
notions we briefly explain how they look in a more classical case of the unit interval. We 
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have that A[o,i] = ^ is the standard Laplacian on [0, 1], and if fJ-[o,i] is the Lebesgue measure 
on [0, 1] then Ar 0il i is self-adjoint and 



£(/,/)= / (/' (x)) 2 dx = - / fA m fdn m 

JO J [0,1] 

for functions in the domain of the Dirichlet or Neumann Laplacian. The energy can also 
be defined by £(/,/) = lim^ 2" E x>y ev n , x ~ y (/(*) - f(v)f where V n = {k/2^ =Q . To 
compute the spectrum of — Am ii one can use the spectral decimation method with inverse 
iterations of the polynomial R{z) = z(4 — z). Each positive eigenvalue can be written as 
A = lim m _^. 0O 4 m A m for a sequence {A m }^ =mo such that A m = R(X m+1 ) and A TOo G {0,4}. 
Then ^R(z) = lim^oo R ok (A~ k z) = 2 — 2cos( v /i) satisfies the functional equation R(y\(z)) = 
5^(42;). In this case a (-A [0 ,i]) C *H _1 {0,4}, the multiplicity is one, and is in the Neumann 
spectrum but not in the Dirichlet spectrum. The eigenfunction extension map is 

i>v,\(%) — cos(\/A \x— v\) - sin(\/A \x— v\) 

sin(v A) 

where v is or 1. 

For much more information on this example and its relation to quantum graphs see [30] 
and references therein. 

2.3. Underlying graph assumptions and Sierpinski fractafolds. Let To = (Vo, Eq) be 
a finite or infinite graph. To define a Sierpinski fractafold, we assume that To is a 4-regular 
graph which is a union of complete graphs of 3 vertices. It can be said that r is a regular 
3-hyper-graph in which every vertex belongs to two hyper-edges. A hyper-edge in this case is 
a complete graphs of 3 vertices, and we call it a cell, or a 0-cell, of r . We denote the discrete 
Laplacian on r by Ar . (In principle, these assumptions can be weakened; see Section [6] 
and Figure EH] for instance). 

Let SG be the usual compact Sierpinski gasket (see Figure |2~2]) . We define a Sierpinski 
fractafold 5 by replacing each cell of T by a copy of SG. These copies we call cells, or 
0-cells, of the Sierpinski fractafold Naturally, the corners of the copies of the Sierpinski 
gasket SG are identified with the vertices of r . 

A fractafold is called infinite if the graph T is infinite. In particular, finite fractafolds are 
compact and infinite fractafolds are not compact. All the details can be found in |37J. In 
this paper we use the same notation as in [37J as much as possible (see also [40J). Since the 
pairwise intersections of the cells of the Sierpinski fractafold $ are finite, we can consider the 
natural measure on the Sierpinski fractafold which we also denote \x. Furthermore, since 
Asg is a local operator, we can define a local Laplacian A on the Sierpinski fractafold as 
explained in 



2.4. Eigenfunction extension map on fractafolds. For any v G Vo and A ^ H ex t there 
is a unique continuous function ip Vi \(-) : 3 — > M such that 

(1) the support of ip v \ is contained in the union of of the cells of the Sierpinski fractafold 
$ that contain v; 

(2) = 1; 

(3) the pointwise eigenfunction equation 

-Aij) Vt x(x) = \ip v ,\(x) 
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holds at any point x G 3\V . 
For any function / on T (and any A as above), we define the eigenfunction extension map 
by 

(2.7) y x f (x) = Y,fo(v)iJ v Ax). 

veVo 

By definition, / = ^xfo is a continuous extension of fo to the Sierpihski fractafold $ which 
is a pointwise solution to the eigenvalue equation above for all x G $\V . Moreover, it is 
known that if fo is a pointwise solution to the eigenfunction equation — Ar /o = Ao/o on r , 
and A ^ {0, 6}, then / = \&a/o is a continuous extension of f to the Sierpihski fractafold # 
which is a pointwise solution to the eigenvalue equation above for all sej. Note that here 
we have A G 9 ! l _1 (Ao), where is as above. The eigenfunction extension map is explained 
in [39] on page 69. 

It is easy to see that ty\ : £ 2 {V ) — > L 2 (^,fi) is a bounded linear operator for any A ^ 
9 t l~ 1 {2, 5, 6}, and its adjoint ty* x : L 2 ($, jj) — > £ 2 {Vo) can be computed as 

(2.8) (ng)(v) = Jg(x)1> ViX (x)dii(x). 

2.5. Spectral decomposition (resolution of the identity). We suppose that the self- 
adjoint discrete Laplacian A ro on T has a spectral decomposition (resolution of the identity) 



(2.9) -A ro = / XdE ro (X). 
which has a form 

(2.10) -A ro /o(«)= / Xj2 p r (\u,v)f (u)dmr (X) 

where m(-) is a spectral measure of —A which is a Borel measure on <r(— A ro ) (see Section [3] 
for more detail). 

We define a function M(A) as the infinite product 

(2.11) M(A) = n niMnini 

where 



A = lim 5 m A m 

m— >oo 



and X m = R(X m+ i). The function M(-) is known from [271 Lemma 2.2 and Corollary 2.4] 
where it appears when the L 2 norm of eigenfunctions on the Sierpihski gasket is computed. 
This function does not depend on the fractafold, but only on the Sierpihski gasket. 
We denote 

= 5 ( 9TH2} U |J 5^~ 1 {3, 5} 

and 

S'oo = 5( |J 5} ) CS C 



m=0 



\m=l 



s 
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Note that for the difference of these two sets we have 

Eoo^ = 3, 5} C DT^O, 6}. 

Theorem 2.3. The Laplacian A is self-adjoint and 

(2.12) 9rV(-Ar )) U c <r(-A) c 9TV(-Ar )) U E M . 
Moreover, the spectral decomposition 

-A = / XdE(X) 

<t(-A) 

can 6e written as 

(2.13) -A= y AAf(A)^d(s ro (*H(A)))* A + £ A£{A}. 

i/ere -E{A} denotes the eigenprojection if X is an eigenvalue (the eigenprojection is non-zero 
if and only if X is an eigenvalue). 

All eigenvalues and eigenf unctions of A can be computed by the spectral decimation method 
as so called offspring of either localized eigenfunctions on approximating graph Laplacians, 
or of eigenfunctions on r . Furthermore, the Laplacian A on the Sierpinski fractafold $ has 
the spectral decomposition of the form 

(2.14) -A/(a:)= J X ^J^P(X,x,y)f(y)dfi(y) S jdm(X) + ^ XE{X}f(x) 
where m = mp o 91 and 

(2.15) P(X,x,y) = M(X) £ ^ v , x (x)^ u , x (y)P ro (m(X),u,v). 

u,v SVb 

Proof. Let T = (Vo, £xi) be as above and let I\ = (Vi, be a graph obtained from T by 
replacing each cell of Tq with the graph shown below. 




The three vertices of the biggest triangle in the above graph replace the three vertices of 
each cell of r . We repeat this procedure recursively to define a sequence of discrete approx- 
imations V n to the fractafold the Sierpinski fractafold On each V n we consider discrete 
energy form, which converge as n — > oo with the same normalization as in Subsection 12.11 
In the limit we obtain a resistance form £ of the Sierpinski fractafold $ and one can use 
the theory of resistance forms of Kigami (see [THJ [19]) to define the weak Laplacian A on 
the Sierpinski fractafold More precisely, the resistance form is a regular Dirichlet form 
on L 2 ($,n) by [T9J Theorem 8.10], for which a self-adjoint Laplacian A is uniquely defined 
(see [19, Proposition 8.11]). One can easily see that in this case the set of continuous com- 
pactly supported functions in DomA such that Af is also continuous (and also compactly 
supported) form a core. For any such function / the Laplacian Af can be approximated by 
discrete Laplacians, that is Af(x) = lim^oo 5 n A n f(x), where A n is the graph Laplacian on 
V n . The limit is pointwise for each x G = (J V n , and is unform on compact subsets of the 
Sierpinski fractafold 5 provided Af is continuous with compact support. The pointwise and 
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uniform convergence of discrete Laplacians in this case is justified in the same as way in the 
case of the Laplacian on the Sierpihski gasket. 

Using notation of Subsections 12.11 and I2.2[ we denote m n = mr o d\ n and 

P n (X,x,y) = M n (X) ^vA^u,x{y)Pr o m\),u,v) 

u,v(zVo 

where M n (A) is defined as the partial product in the definition of M(A). We further denote 

/ n-1 ^ 

E n = 5 hH^{2}U |J 5m 9C{3,5} 

V m=0 / 

and let E n X be the eigenprojection of — A n corresponding to A. Then we have the discrete 
version of the formula (I2.14p because of the computation in P, Theorem 3.3] (see also Sections 
[3] and H] below, where Pr (X,u,v) is denoted by P\(u,v)). Note that in [H Theorem 3.3] the 
normalization factor was l/(0i?'), where 4>(z) = ^_^i_ 2z ) an d R(z) = z(5 — 4z), which 
produces the normalization factor 

(5 - 4z)(l - 2z) _ 1 (1 - 4z/5)(l - Az/2) 

(3-2;z)(5-8z) ~ 3(1 -4z/6)(l -8z/5)' 

which is the same as in (12. lip . Here 4z replaces A m because the distinction between proba- 
bilistic and graph Laplacians, and the extra factor | is because of the integration in (12. 14ft . 

Let u and / be continuous functions on the Sierpihski fractafold $ with compact support 
and let 

v = (-A + l)- 1 f. 

The usual energy and L 2 estimates imply that v G Dom(A) is continuous, square integrable, 
and —Av = f — v. We have, by the discrete approximations, that the inner product (u, v)i? 
is equal to 



£R-i(<x(-Ar ))\S o 

and so we have the relation 



(u, v) L 2= I (u, dE(X)f) L 2 

Ja(-A) A + 1 

when u, f are continuous functions with compact support. The theorem then follows by the 
general theory of self-adjoint operators [32} Section VIII. 7]. □ 

2.6. Infinite Sierpinski gaskets. As a collection of first examples we consider the infinite 
Sierpihski gaskets, where the spectrum was analyzed in [3llIIJ[3T]- First, note that up to a 
natural isometry there is one infinite Sierpihski gasket with a distinguished boundary point 
(and hence it is not a fractafold), and there are uncountably many non-isometric infinite 
Sierpihski gaskets which are fractafolds (see [H] for more detail). 

If an infinite Sierpihski gasket fractafold is build in a self-similar way, as described in 
[36| l4Tj . then the spectrum on r is pure point with two infinite series of eigenvalues of 
infinite multiplicity. One series of eigenvalues consists of isolated points which accumulate 
to the Julia set 3r of the polynomial R, and the other series of points are located on the 
edges of the gaps of this Julia set (the Julia set in this case is a real Cantor set of one 
dimensional Lebesgue measure zero). The set of eigenvalues E on r consists of 6 and all 
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Figure 2.3. An illustration to the computation of the spectrum on the infi- 
nite Sierpihski gasket. The curved lines show the graph of the function JH(-), 
the vertical axis contains the spectrum of a(— Ar ) and the horizontal axis 
contains the spectrum cr(— A). 



the preimages of 5 and 3 under the inverse iterations of R. In this case formula (I2.14p is 
the same as the formulas for eigenprojections in [41] . The illustration to the computation 
of the spectrum in Theorem 12.31 is shown in Figure 12.31 where the graph of the function 

is shown schematically and the location of eigenvalues is denoted by small crosses. The 
spectrum cr(— A) is shown on the horizontal axis and the set of eigenvalues So of — Ar is 
shown on the vertical axis. 

A different infinite Sierpihski gasket fractafold can be constructed using two copies of an 
infinite Sierpihski gasket with a boundary point, and joining these copies at the boundary. 
This fractal first was considered in [2], and has a natural axis of symmetry between left and 
right copies. Therefore we can consider symmetric and anti-symmetric functions with respect 
to these symmetries. It was proved in [UJ that the spectrum of the Laplacian restricted to 
the symmetric part is pure point with a complete set of eigenfunctions with compact support. 
For the anti-symmetric part the compactly supported eigenfunctions are not complete, and 
it was proved in [31] that the Laplacian on To has a singularly continuous component in 
the spectrum, supported on 3r, of spectral multiplicity one. As a corollary of these and our 
results we have the following proposition. 
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Proposition 2.4. On the Barlow-Perkins infinite Sierpinski fractafold the spectrum of the 
Laplacian consists of a dense set of eigenvalues £H _1 (£o) of infinite multiplicity and of a 
singularly continuous component of spectral multiplicity one supported on £H -1 (£[_r). 

3. General infinite fractafolds and the main results 

Consider a fractafold with cell graph T, so T is an arbitrary infinite 3-regular graph. The 
spectrum of — Ar is contained in [0,6], and by the spectral theorem there exist projection 
operators Ej corresponding to intervals / C [0,6]. Because we are in a discrete setting we 
can say a lot more. There is a kernel function Ej on V x V such that 

(3.1) E I f{a) = Y,Ei{a ) b)f{h) 

bar 

and I — > Ej(a, b) is a signed measure for each fixed a, b. Since there are a countable number 
of such measures, we can find a single positive measure /i on [0,6] such that 

(3.2) £/M) = Jp x (a,b)dfi(X) 

for a function P\(a,b) defined almost anywhere with respect to fi (so P\(a,b) is just the 
Radon-Nykodim derivative of Ej(a, b) with respect to /i). In fact, by a theorem of Besicovitch 

3.3 P A (a, b) = hm — l —— 

<^o n([\ - e, A + ej) 

for [M — a.e.X. (if fi is absolutely continuous this is just the Lebesgue differential of the integral 
theorem). It follows from f)3.3p that 

(3.4) -A r P A (- s 6) = AP A (.,6) 
for fi — a.e.X. Thus if we define the pointwise projections 

(3.5) P A /(a) = ^P A (a, b)f(b) 

then the spectral resolution is 

(3.6) / = J^Pxfdfi(X) 
with 

(3.7) - A r P x f = XP x f, 

where E C [0, 6] is the spectrum. In other words, ( 13. 6 p represents a general function / 
(we may take / G ^ 2 (r), or more restrictively a function of finite support) as an integral 
of A-eigenfunctions. Note that typically P A / is not in £ 2 (T). Also, the measure \i and the 
kernel Pa are not unique since one may be multiplied by g(X) and the other by for any 
positive function. We are not aware of any way to make a "canonical" choice to eliminate 
this ambiguity. 

We also observe that the measure /i does not have a discrete atom at A = 6. In other words, 
there are no i 2 {T) 6-eigenfunctions. Indeed, for a 3-regular graph, there exist 6-eigenfunctions 
if an only if the graph is bipartite, in which case the 6-eigenfunction alternates ±1 on the 
two parts. Since we are assuming T is infinite, this eigenfunction is not in £ 2 (T). 
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Let T denote the edge graph of T. Then T is 4-regular. Let A ro denote its Laplacian. 
Define 

(3-8) P x (x,y) = -^-J2J2 P ^V 

(there are 4 terms in the sum). Let Eq denote the space of 6-eigenfunctions in £ 2 {Tq) (this 
may be 0) and write Pq for the orthogonal projection of £ 2 (To) onto Eq. 

Theorem 3.1. The spectral resolution of — Ar is given by 
(3.9) F = P 6 F + [ P x Fd/i(X) 



where 

(3.10) - A To P x F = XP X F 
for \i — a.e.X, and 

(3.11) P x F(x) = J2 p ^y) F (y)- 

yer 

In particular, spect(— Ar ) = X (if E 6 = 0) or S U {6}. 

For the proof we require some lemmas. 
Following jlQ] we define the sum operators 

Si : f(T ) 

and 
by 

(3.12) SJ(x) = f(a) + fib) if x is the edge (a, b) 
and 

(3.13) 5 2 F(a) = F(x) + F(y) + F(z) if x, y, z are the edges containing a. 

Lemma 3.2. S2S1 = 6/ + Ar and S1S2 = 61 + Ar . in particular, S2S1 is invertible, Si is 
one-to-one and S% is onto. 

Proof. The formulas for S2S1 and S1S2 are simple computations. Since there are no 6- 
eigenfunctions in £ 2 {T), we obtain the invertability of S2S1 (see also |1Q]). □ 

It follows from Lemma O that E 6 = (ImSi) 1 - and £ 2 {T ) = ImS'i © ^e- 

Lemma 3.3. For any A / 6, — A-pf = Xf if and only if —A^Sif = XSif. In particular, 
sp(-A ro ) = sp(-A r )U{6}. 

Proof. Suppose — A r / = Xf. Since — A r = 61 — S 2 Si we have — S^Si/ = (A — 6)/. Apply 
Si to this identity and use — Ar = 61 — S1S2 to obtain — Ar Sif = XSif. Similarly, we can 
reverse the implications. Note that the condition A 7^ 6 implies that Sif is not identically 
zero (see also |1Q]). □ 
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Lemma 3.4. Let F G £ 2 {^o) be orthogonal to Eq (if Eq is nontrivial) . Then F = Sif for 

(3.14) f = (QI + A r )- 1 S 2 F 1 . 

Moreover we have 

(3.15) P x = -L-SxP x S 2 . 

o — A 

Proof. For / defined by (13.41) we have S 2 Sif = S 2 F by Lemma [3T2l Since S 2 is injective on 
Eq and Sif G Eq we conclude Sif = F. 

By definition, P x F(x) = 7 -Fa (a, b)F(y) and this is equivalent to (13 .150 by 

yeTo a£x bey 

the definition of Si and 5*2. □ 

Proof of Theorem \3.1[ It suffices to establish (13.91) for F G Eq. For / defined by (13.141) . we 
apply Si to (13.61) to obtain 

F = J S 1 P x fd f i(\) = J S 1 P x (6I + A G )- 1 S 2 Fd f i = j -L-S 1 P x S 2 Fd t M(X) 

since P A (6J + A r )- X = ^P A . Then (E3]) follows by ( 13~T5|) . We obtain (13401) from (13T7D 
and Lemma [3.31 □ 

In order to give an explicit form of the spectral resolution for any particular T, we need 
to solve two problems: 

(a) Find an explicit formula for P\(a, b); 

(b) Give an explicit description of Eq and the projection operator Pq. 

In addition, there is one more problem we would like to solve in order to obtain an explicit 
Plancherel formula. We can always write 

(346) ||/||| 2(r) = J < P x fJ > dn(\) 

and 



(3.17) I |F| | £ 2 2(ro) = \\PeF\\2 + j^<PxF 1 F>dfi(X) 

for a reasonable dense space of functions / and F (certainly finitely supported functions 
will do). What we would like is to replace < P x f, f > and < P X F, F > by expressions only 
involving P x f and P X F and some inner product on a space of A-eigenfunctions. Note that 
from (I3.2p and the fact that Ei is a projection operator we have 

(3.18) < P x f, f >= \ima([X - e, A + e])|| \ E [x ^ x+e] f\\ 2 2 

<^o /i{[A — e, A + e\) 

for /j, — a.e.X. This suggests the following conjecture, 

Conjecture 3.5. For /i — a.e.X there exists a Hilbert space of X-eigenf unctions £ x with inner 
product <, > A such that P x f G £a f or H — a.e.X for every f G i 2 (T), and 

(3.19) <P x f,f>=<P x f,P x f> x . 

Moreover a similar statement holds for < P X F, F > . 

Our last problem is then 
(c) Find an explicit description of £ x and its inner product, and transfer this to £ x of r . 



14 



ROBERT S. STRICHARTZ AND ALEXANDER TEPLYAEV 



4. The Tree Fractafold 

In this section we study in detail the spectrum of the Laplacian on the tree fractafold 
TSG (Figure 14.11 ) whose cell graph T is the 3-regular tree. In a sense this example is the 



Figure 4.1. A part of the infinite Sierpihski fractafold based on the binary tree. 

"universal covering space" of all the other examples, if we "fill in" all copies of SG with 
triangles. 

We begin by solving problem (b). 

Lemma 4.1. For any fixed z in T define 

F z (x) = ^-\) d{x ' z) ■ 
Then F z e £ 2 (T ) with \ \F z \\p {To) = 1 and F z e E 6 . 

Proof. Note that z has 4 neighbors {y%, y^ 2/3, 2/4} in r with d(yj, z) — 1, so 

4 

-A To F z (z) = 4F M (z)-Y,HVi) 

verifying the 6-eigenvalue equation at z. 

On the other hand, if x 7^ z then the 4 neighbors {yi, y^-, 2/3, 1/4} of x may be permuted so 
that d(yi, z) = d(x, z) — 1, d(y 2 , z) = d(x, z), and d(y 3 , z) = d{y±, z) = d(y, z) + 1. It follows 
that 

4 1 
-A ro F 2 (x) = 4F z (x) - ^ F z {y 3 ) = F z (x)(A - (-2 + 1 - 2 • -)) = 6F z (x) 

verifying the 6-eigenvalue equation at x. Finally 

ii^ii,V«) = ^ 1 + 4 -4) 2 + 8 -(i) 2 + ---) = ^(i + i + ^ + ^ + ---) = i 

(See Figure S3). □ 
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Figure 4.2. Values of y/3F z (the center point is z). 

Remark 4.2. It is easy to see from the 6-eigenvalue equation that F z is the unique (up to a 
constant multiple) function in E Q that is radial about z (a function of d(x, z) ). 

Lemma 4.3. ^F z (x)F y (x) = V3F z {y). 

X 

Proof. Fix z. Then the left side is a 6-eigenfunction of y and is radial about z, so it must be 
a constant multiple of F z (y). To compute the constant set y — z, and the left side is 1 while 

f z {z) = -fe. a 

Definition 4.4. Let P e (x,y) = ^F x (y) = |(— ^ d ( x < y > and define the operator 

(4.1) P 6 F(x) = Y,Mx,yMy). 

y 

Theorem 4.5. P§ is the orthogonal projection ^ 2 (Fo) — > E§. 

Proof. Lemma [4.31 shows PqF z = F z . Now we claim that the functions F z span Eq. Indeed, 
if F is in E Q and is orthogonal to F z , then we can radialize F about z to obtain a function F 
that is still in E 6 and orthogonal to F z . Since F must be a multiple of F z it follows that it is 
identically zero. Since F(z) = F(z) it follows that F(z) = 0. Since this holds for every z, we 
have shown that the orthogonal complement of the span of F z is zero. This shows Pq is the 
identity on E$. Also PqEq = by the orthogonality of different parts of the spectrum. □ 

Note that {F z } is not an orthonormal basis of Eq, since < F z , F y >= \/3F z (y) by Lemma 
14.31 The next result shows that it is a tight frame. 

Theorem 4.6. For any F 6 Eq we have 

(4-2) £|<F,F,>| 2 = 3||F||, 2 2(ro) 

Z 

Proof. We may write F = ^a{y)F y . Then | | 2 2(ro) = S ^ S ^a{y)a(z)\/?>F z (y). But 

y y z 

< F, F z >= a{y)V3F z {y) and so 
y 

J2\<F,F Z >\ 2 = 3j2J2J2a(y)'aJ^)F y (z)F y/ (z) 

Z z y y' 

(4-3) = 3YYa(y)aJ^F yl (y) =3\\F\\ 2 e2(ri)) 

y y' 
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□ 

It follows from polarizing ( 14. 3 p that we may also write PqF = | < F, F z > F z . 

z 

The solution of problem (a) is due to Cartier [4]. We outline the solution following [9]. 

Definition 4.7. Let z e C with 2 2z ~ l ^ 1. Let c(z) = ~ , c{\ - z) = 1 2 -/S 2 ^-i and 
<p z (n) = c(z)2~ nz + c(l - ^)2- n ( 1 - 2 ). 

Remark 4.8. Note thatc(z) and c(l—z) are characterized by the identities c(z)+c(l — z) = 1 
andc{z)2~ z +c(l-z)2 z - 1 = c(z)2 z + c(l- z)2^ z which imply <p z (0) = 1 andcp z {l) =<p z (-l). 

Theorem 4.9. For any fixed y G T ; let f y {x) = tp z (d(x,y)). Then 

(4-4) -A r /, = (3-2^-2 1 - z )/, 

and f y may be characterized as the unique (3 — 2 Z — 2 1 ~ z )-eigenf unction that is radial about 
y and satisfying f y (y) = 1. 

Proof. Uniqueness follows from the eigenvalue equation. To verify the eigenvalue equation 
we do the computation separately for x ^ y and x = y. For x ^ y note that x has 
two neighbors, x\ and x%, with d(xi,y) = d(x2,y) = d(x,y) + 1 and one neighbor, X3, 
with d(x-s,y) = d(x,y) — 1 so the eigenvalue equation is immediate. On the other hand y 
has three neighbors, Xi,x 2 ,x 3 , with d(xj,y) = 1, and the eigenvalue equation follows from 
ip z {l) = ip z {-\). " □ 

Note that there is no choice of z that will make f y belong to £ 2 {T). However, the choice 
z = I + it gets close. Indeed \ipi +it (d(x, y))\ 2 ~ 2 n ■ 2~ n just diverges. So it is natural 

n 

to conjecture that these eigenfunctions give the spectral resolution of — Ar on £ 2 (T). In fact 
the following proposition is the content of Theorem 6.4 on p. 61 of [9]. 

Proposition 4.10. By periodicity we may restrict t to < t < Write X(t) = 3 — 

2 v / 2cos(tlog2) = 3-2^ +i *-2H* andJ2 = [3-2^,3 + 2^] « [0.17,5.83] C [0,6]. Define 

(4-5) Ptf{x)=Y,<Pi + uMx,y))f{y). 

y 

Note that -A r P t f = X(t)P t (f). Then 

(4.6) f(x) = / l0S2 PJ(x)dm(t) 

Jo 

for the measure 



(4.7) dm(t) = ^ 

671 



C(- + it) 



^ (My; 1 

tt(1 + 2 sin 2 (t log 2)) 



It is convenient to change notation so that the eigenvalue A rather than t is the parameter, 
e easily co: 
Note that 



We easily compute t = cos 1 (§71) ■ 



d\ = 2v / 21og2sin(tlog2)c/t, sin 2 (tlog2) = (-A 2 + 6A - l)/8, 

and 

1 + 2 sin 2 (t log 2) = (-A 2 + 6A + 3)/4. 



SPECTRAL ANALYSIS ON INFINITE SIERPINSKI FRACTAFOLDS 17 



If we write Pa — Pt then the spectral resolution is 

/•3+2 v / 2 

f(x) = / P x f(x)dm(X) 



3-2V2 



for 



dm(\) = —= — d\. 

V2n(-\ 2 + 6A + 3) 



Now suppose F e ^ 2 (r ) lies in E$. Then we may write F = Sif for / = (6/ + Ar)^ 1 S 2 F 
in £ 2 (T). Indeed we know that 6 is in the resolvent of — Ar so / is well-defined, and then 
S2S\f = S2F by Lemma I3.2L Since S 2 is injective on E§ and S\f G Eq we conclude 
S 1 f = F. 

By Proposition 14.101 we have 

(4.8) Stf = f SxPxfdmiX), 



and of course — Ar SiP x f = XSiP x f by Lemma [3731 so we define P\F = SiP\f and we 
obtain the spectral resolution of F: 

(4.9) F = J P x Fdm(\). 

Note that P x {61 + Ar)" 1 = ^P x so P X F = ^StP^F. 
We may write this quite explicitly as follows: 

Lemma 4.11. Define 

(4.10) Vz(n) = c{z)2~ nz + 5(1 - z)2~ n ^ 

for c(z) = (2 + 2~ 2 + 2 z )c{z). Note that ip z (n) = 2y z (n) + cp z (n + 1) + (f z (n - 1). Then 

(4.11) S x P x S 2 F{x) = l^i+uidfayVFiv). 

y 

Proof. S 2 F(b) = J2 F (y">- There 

are 3 terms in the sum, and y ~ b means the edge y joins 
b and one of its neighbors in T. Then we compute 

(4.12) P X S 2 F(a) = X) E ^^t(d(a, b))F(y) 



2 

6er j/~6 



and 

(4.13) S 1 P x S 2 F(x) = ^EE^^^W 

fegr 

where a ~ x means that a is one of the vertices in the edge x. Suppose x 7^ y and let 
n = d(x, y) with n > 1, (Figure S3] shows the To graph for n = 2). 

Then x ~ ai and a; ~ a 2 while y ~ 61 and y ^ b 2 with rf(ai,6 2 ) = <i(a2,&i) = n, 
d(ai, b\) = n — 1, and d(ei2, 62) = n + 1- The result follows in this case. When x = y we have 
d(x,y) = and ai = 61, a 2 = &2 so d(ai,b 2 ) = d(a 2 ,bi) = 1 and d(ai,bi) = d(a 2 ,b 2 ) = 0. 
The result follows because (fi +it (— 1) = (pi +it (l). □ 
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Theorem 4.12. For any F e £ 2 (r ) we have the explicit spectral resolution 

(4.14) F = P 6 F + J P x Fdm{\) 

for 

(4-15) PxF(x) = ^+iMx, y))F(y). 

The Theorem follows by combining Lemma 14.111 and Proposition 14.101 We note that the 
proof of Proposition 14.101 involves an explicit computation of the resolvent (A/ + Ap) _1 for A 
outside the spectrum of — Ar, followed by a contour integral to obtain the spectral resolution 
from the resolvent. We sketch some of these ideas and then show how to carry out a similar 
proof of Theorem 14.121 

On £ 2 {T) we define 

(4.16) H z f{a) = Y J ^ zd{aM f{h). 

b 

A direct computation shows 

(4.17) (\I + A r )H z f = (2- z -¥)f 

for A = 3 - 2 Z - 2 • 2~ z . 

Note that = cosh((,2 — |)log2), and in order to have H z bounded on £ 2 (T) we need 

^tz > ~. This shows spect(— Ar) = £ and (XI + Ar) -1 = 2 _ 2 1 _ 2Z if z for z ^ E. 
On £ 2 (T ) we define 

(4.18) H z F{x) = J2 2~ zd(x ' y) F(y). 

y 

Lemma 4.13. spect(-A ro ) _1 = SU{6} and (A/+A) -1 = 2 . 2 - z 1 _ 2 ,_ l H z forz spect(-A To ). 

Proof. Note that H z is bounded on £ 2 (T ) for ^R,z > |. Also A = 6 corresponds to z = 1 + 
for which 2 • 2~ z — 2 Z — 1 = 2(— |) — (2) — 1 = 0. Now fix x and consider its four neighbors, 
Xi,X2,x 3 ,X4 (so d(x,Xj) = 1). For any fixed y ^ x we may order them so that d(xi,y) = 
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d(x 2 ,y) = d(x,y) + 1, d(x 3 ,y) = d(x,y), d{x^y) = d(x,y) — 1. It follows that 
(XI + A To )H z F(x) = (A - A)H z F(x) + ^ H z F{x 3 ) 

3 

= (A - 4)F(x) + ^ 2- 2 F(s) + (A - 4) ^ 2- zd ^ y) F{y) + ^ ^ 2- zd{x > y) F(y) 
j vfi* 3 y^x 

(4.19) = (2 • 2~ z - 2* - l)F{x) 

and the result follows. □ 
For / G ^ 2 (r), we have 

(4.20) /=-L /(AJ + Ar)- 1 /^ 

27tz y 7 

for any contour 7 that circles £ once in the counterclockwise direction. We choose 7 as 
shown and take the limit as 5 — > + . The contribution from the vertical segments goes to 
zero so 






. 5 . 


5 


£ 


7 

( 









Figure 4.4. The contour 7 for integration in (I4.20p . 

(4.21) f= lim — / ((A-^ + A r )" 1 -(A + ^ + A r )" 1 )/rfA. 

If z = \ + e + it for e > then |=| = cos(t log 2 - ze log 2) and 

(4.22) A = 3-2v / 2cos(tlog2)cosh(elog2) - «2v / 2sinh(elog2) sin(t log2). 

For t > we have A ~ 3 — 2y / 2cos(tlog2) — while for t < we have A ~ 3 — 

2\/2cos(tlog2) + iS with 5 > 0. Thus 

(4.23) 

« A - * + Ar )" " ( A + 15 + Ar )") / = 2 -i-^ 2 i + , g i + ^- 2 -H^- 2 i-^ -^ 

so we obtain 
(4.24) 



^ 2vri 



log 2 



T-7 H h+itf 



2~2- u - 2? +it 2 



— Hi_ it f) 2v / 21og2sin(tlog2)dt. 



This is the same as / = J log2 Ptfdm(t). 
For F G ^ 2 (r ) we have 



(4.25) 



F=^- /(A/ + ArJ- 1 ™ + / (A/ + ArJ- 1 ™ 
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where 7 is as before and 7' is a small circle about 6. Taking the limit we obtain 
F = lim -L / ((X-iS + ArJ- 1 ^- (X + iS + A ro )- 1 F)d\ 

(4.26) + lim / (6 + 5e ie + Ar^Fide^dO. 

8^0+ 2m J 

As before we can write the first term as 

v / 21og2 r ( 1 

tt« is \25 _it - 2t+ rt - 1 2 2t +it - 2t" rt - 1 



( 4 - 27 ) JL — — / ( ^ _ 9 ^ t _ : " l^r—^T—^_ lt F ) sin(t log 2)dt. 



which we identify with J s P\Fdm(X), while the second term is P%F. 

Next we discuss an explicit Plancherel formula on T, given in terms of the modified mean 
inner product 

(4.28) <f,g> M = Jim^ 

d(x,xo)<N 

We will deal with eigenspaces for which the limit exists and is independent of the point xq. 
Note that this is not the usual mean on T, since the cardinality of the ball {x : d(x, xo) < N} 
is 0(2 n ), but it is tailor made for functions of growth rate 0(2 -d ( x ' x °)/ 2 ), which is exactly 
the growth rate of our eigenfunctions. 

We expect that analogous results are valid for A;-regular trees for all k, but to keep the 
discussion simple we only deal with the case k — 3 that we need for our applications. 

Lemma 4.14. For all n and t 

(4.23) ,, +(t(ll) ^(3e„ s(lltl o g2 ) + «|) 2 - 

Proof. From the definition, 

n+it (n)=^(c(l + it)2-^ 2-/ 2 . 

The result follows from the explicit formula for c(^+it) and some trigonometric identities. □ 

In what follows we write ip for V?±+it to simplify the notation. 
Lemma 4.15. Let 

Then for any integers k and j , 

N N 

(4.31) lim -Y,2 n+ h(nMn + k) = lim - £ 2 n+ ^ <p(n + J + k)<p(n + j) 



jV->oo N N^oo N 

n=l n=l 



jU(A) cos(A;t log 2). 
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Proof. It is easy to see that (I4.3ip is independent of j, so we take j = 0. Then by (I4.29p 

2 n+ ^(n)<p(n+k) = - (3 cos(nt log 2 ) + sm ( nt log S ) (3cos(ntlog2) cos(fctlog2) 
9 \ tan(tlogz) / 

sin(ntlog2)cos(fctlog2) cos(ntlog2) sin(fet log 2) ' 

— 3sin(ratlog2) smlkt log 2) H — - — 

v 6 ; v B ' tan(tlog2) tan(tlog2) 

Now use the following identities 

1 N 1 N 1 

lim — > cos 2 na = lim — > sin 2 na = - 

N^oo N ^ N^oo N Z — ' 2 

n=l n=l 

and 

1 N 

lim — > cos na sin na = 

7V->oo N ^— ' 
n=l 

to see that the limit in fl4.31f) equals 

1 / /7 , „, 3 sinolog 2) 3 sinolog 2) cos(fctlog2)\ 1,,,. fl , „, 

— 9 cos (Art log 2 + ; - 6 / ) - V + — v 9/ - 6 { = — 6(A) cos(kt log 2). 

18 V tan(tlog2) tan(tlog2) tan 2 (tlog2)y 18 v ; v B J 

□ 

Lemma 4.16. For any A in the interior ofH and x\ G V, < P\5 X1 , P x 6~ Xl >m exists and is 
independent of the base point x , and 

(4.32) < P x 5 Xl ,P x 5 Xl > M = ^b(X). 

Proof. P x 5 Xl (x) = <p(d(x, xi)) and <p(n) = 0(2~ n / 2 ) by (jQ9l . It follows easily that the limit, 
if it exists, is independent of the choice of x , since if d(x , x' Q ) = k then B n _ k (x' Q ) C B n (x ) C 
fi n+ fc(x' ), and the division by N in (14.281) makes the difference go to zero as iV — > oo. We 
will prove the existence of the limit by computing (14 . 3 2 [) with Xq — X\ . 

Note that there are exactly 3 • 2 n_1 points x with d(x, Xq) = n for n > 1, and we can ignore 
the point x = x\ in computing the limit. Thus 

3 " , 1 



< P X 5 X11 P X 5 X1 > M = lim _^2>(n) 2 = -6(A) 

V-S.OO ZiV z — ' 12 
n=l 

by Lemma 14.151 □ 



Lemma 4.17. Suppose d{x\,X2) = k and A is in the interior of E. TTien < P X S X1 , P\8 X2 >m 
exists and is independent of the base point xq, and 

(4.33) < P X 5 X1 ,P X 5 X2 > M = ^b(X)(p(k). 

Proof. The proof of independence of the base point is the same as in Lemma I4.16[ so we 
compute the limit for xq = x\. Except for a few points when n is small that don't enter into 
the limit, we may partition the points with d(x,xi) = n as follows: 
2™ points with d(x, x 2 ) = n + k, 

2™-i- 1 points with d(x, x 2 ) — n + k — 2j for 1 < j < k — 1, 
2 n ~ k points with d(x, x k ) = n — k. 
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d(x, £2) = n + k — 2 d(x, Xq) = n — k + 2 



d(x, x 2 ) = n + k 





d(x, x 2 ) = n 



Figure 4.5. Partition of points x with d{x,x\) = n. 



This implies 

N 



N / k-l 

n 53 ( 2 MnMn + k) + -J2 2""V(«)^(n + k - 2j) + 2 n ~ k (p(n)<p{n - k) 



71=1 



3=1 
fe-1 



= ^6(A)2~ fc / 2 ^cos(H log 2) + 5 53 cos ^ ~ 2 ^')* lo S 2 + cos ( kt lo S 2 ) 
by Lemma 14.151 



k-l 



However, the trigonometric identity sin(a) 53 cos (^ ~ 2j)a = sin(fca) cos(a) implies 

i=o 
fe-i 

2 cos(A;t log 2) + § 53 cos ( fc ~~ ^ lo S 2 
3=1 
jfe-i 

= § cos(fct log 2) + § 53 cos ( fc ~ 2 ^)* lo S 2 



j=o 

= § (cos(Hlog2) + =f^)2^ 
by Lemma I4.14| which implies (I4.33P . 
Theorem 4.18. Suppose f has finite support. Then 

(4.34) < P x f, f >= 126(A)- 1 < P x f, P x f > M . 
Proof. Since < P\S Xl ,5 X2 >= ip(d(xi, x 2 )) we can rewrite f |4.33|) as 

< P\5 X1 ,5 X2 >= 126(A)" 1 < P x 5 Xl ,P x 5 Xl > M 
and (I4.34p follows by linearity. 

Corollary 4.19. For f G £ 2 (T), for /i a.e. X, < P x f,P x f >m exists, and 

(4.35) ||/|| 2 2{r) = J < Pxf, Pxf > M 126(A)- 1 rf / x(A). 
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Proof. For / of finite support, (I4.35P follows from (I4.34p and (I3.16p . It then follows for 
/ G £ 2 (T) by routine limiting arguments. □ 

To complete the solution of problem (c) for this example we need to transfer the result 
from T to r . Define the modified mean inner product on T by ( 14 . 2 8 [) again, where / and 
g are functions on T and x and Xq vary in T . 

Lemma 4.20. For any integers k and j, 

1 N k 1 N k 

lim —\^2 n+ ^ip(n)if)(n + k) = lim — 2 n+j+ ^ip(n + j)ip{n + j + k) 

n=l n=l 

(4.36) = ^z|)i&(A)cos(A;tlog2). 

Proof. As in the proof of Lemma [4. 151 it is clear that (I4.36j) is independent of j, so we may 
take j = 0. Since ip{k) = 2ip{k) + ip{n — 1) + ip{n + 1) we may reduce (I4.36P to (I4.3ip as 
follows: 

N 

lim i 2 n+ ^ip(n)tp(n + k) 

n=l 

N 

k 



lim i V2 n+ 2 (2ip(n) + <p(n-l) + <p(n+l)) {2<p{n+k) + ip{n+k-l) + ip{n+k+l)) 

n=\ 

f [(4 + 2 + |) cos(H log 2) + 2( ^2 + j= 



n=l 



[(4 + 2 + |) cos(fct log 2) + 2( y/2 + 4=) (log ( A; + l)i log 2 + log(fc - l)t log 2) 



+ cos(A; + 2)t log 2 + cos(A; - 2)t log 2] 
^cosfctlog2[(4 + 2 + |) +^( v / 2 + ^)costlog2 + 2cos2tlog2] 
= ^ cos kt log 2(-| + 2 cost log 2f 
and (jQSD follows since ^~ + 2 cos t log 2 = ^1. □ 

Lemma 4.21. For an?/ A in the interior o/E and xi G r 0; < P\d~ Xl , P\d~ Xl >m exists and is 
independent of the base point x , and 

(4.37) <P X 6 X1 ,P X 6 X1 >M= h -^r- 

162 

Proof. The proof that the limit is independent of the base point is the same as in Lemma f4. 161 
so we compute (I4.36j) with xq = X\. Note that for n > 1 there are exactly 4-2 n_1 points x in Vq 
withd(x,xi) — n. For such points P\5 X1 (x) = gr^|^(n) = fzt| (2ip(n) + ip(n— 1) + ip(n+l)) 
and so 

1 2 1 °° 

< P x 5 Xl ,P x 5 Xl > M = —-— ■ - Jim - 2™ (2^(n) + ^(n - 1) + <p(n + l)) 2 

^ ' n=l 

and (jQ7|l follows from (OB) . □ 

Lemma 4.22. Suppose d{x\,Xo) = k and A m t/ie interior o/E. T/jen < P X 5 X1 , P X 5 X2 >m 
exists and is independent of the base point, and 

(4.38) < P X S X1 ,P X S X2 > M = ^ ■ 3 J_ x M k)- 
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Proof. As before we can take the base point Xq = x\. For n > k we can sort the 2 n+1 points 
x with d(x, Xx) = n as follows: 
2 n points with d(x, X2) = n + k, 

2 n ~i points with d(x, X2) = n + k — 2j + 1 for 1 < j < k, and 
2 n ~ k points with d(x, X2) = n — k. 



d(x, x 2 ) = n + k — 3 



d(x,x 2 ) = 




d(x, x 2 ) = n + k — 1 d(x, x 2 ) = n — k + 1 



Figure 4.6. Partition of points x with d(x,x\) = n. 



Thus we have 



11 1 N 

< PxS Xl , P X 6 X2 > M = ■ - lim - 



(6 - A) 2 9 iv^oo N ^ 

v ' n=l 

i/)(n) ( 2 n tj){n + k) + Y^ 2 n ~ j ^(n + k - 2j + 1) + 2 n "^(n - k) 



9-36 



cos(/rf log 2 ) + "Tf cos ^ ~ 2j + 1 )' log 2 + cos (^ lo S 2 ) 



by (ICTjl . 

To complete the proof we need to show 



2 cos(H log 2) + ^ cos(A; - 2j + l)t log 2 
= 3^(2^(fc) + tp(k - 1) + ^ + 1)). 
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As we saw in the proof of Lemma EH\ <p(k) = §2- fc / 2 (2cos(£;tlog2) + ± ^cos(£;-2j)tlog2) 

3=1 



so 



fc-1 



2<p(k) +p{k - 1) + <p(k + 1) = |2" fc / 2 4 cos(fct log 2 + ^ cos(A; - 2j) log 2 



fc-2 



□ 



+2 v / 2cos(A; - l)Mog2+ ^ ^ cos(/c - 2j - l)Hog2 + v / 2cos(/c + l)tlog2 

j=i 
h 

3=1 

and the result follows by standard trigonometric identities. 
Theorem 4.23. Suppose F has finite support on T . Then 
(4.39) < P\F, F >= 366(A)- 1 < P\F, P X F > M . 

Proof. Since < P\S Xl ,5 X2 >= 3 ,^_ x ^ (d(xi,X2)) we can rewrite (14. 3 7ft as < P\5 Xl ,5 X2 >= 
366(A)- 1 < P\5 Xl ,P x 5 Xl > M and (jCTJ) follows by linearity. □ 

Corollary 4.24. For F 6 £ 2 (T ), for fi-a.e. A ire E, < P A P, P A P > M exzsfo, and 



|2 

l^ 2 (r ) 



6^ ||2 



<P A F,P A F> M 366(A)- 1 rf/i(A). 



Proof. Same as for Corollary 14.191 



□ 



We end this section with a description of 5-series eigenf unctions on the graph Ti (note 
there are no 5-eigenfunctions on the graph r ). One can easily see that on I\ there are no 
finitely supported 5-eigenfunction, there are no radially symmetric 5-eigenfunctions, and that 
5-eigenfunctions do not correspond to cycles. By by an argument similar to Theorem 14 . 5 1 one 




Figure 4.7. A part of Ti with a 5-eigenfunction (values not shown are equal 
to zero). 



can show that eigenfunctions in Figure I4T7I (with their translations, rotations and reflections), 
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are complete in the eigenspace E 5 on Ti. We do not give an explicit formula for the 5- 
eigenprojections on T n . One can see that for each n > 1 there are eigenfunctions on T n that 
resemble those in Figure W7f\ and also finitely supported 5-eigenfunctions (see Remark 15. ip . 

5. Periodic Fractafolds 

Remark 5.1. Note that on a periodic graph, linear combinations of compactly supported 
eigenfunctions are dense in an eigenspace (see [231 Theorem 8], [22J and [2H Lemma 3.5] j. 

The computation of compactly supported 5- and 6- series eigenfunctions is discussed in 
detail in J37J HI], and the eigenfunctions with compact support are complete in the corre- 
sponding eigenspaces. 

In particular, [371 HI] show that any 6-series finitely supported eigenfunction on T n+ i is 
the continuation of any finitely supported function on T n , and the corresponding continu- 
ous eigenfunction on the Sierpinski fractafold # can be computed using the eigenfunction 
extension map on fractafolds (see Subsection \2.J$ . Similarly, any 5- series finitely supported 



eigenfunction on T n+1 can be described by a cycle of triangles (homology) in Y n , and the 
corresponding continuous eigenfunction on the Sierpinski fractafold $ is computed using the 
eigenfunction extension map on fractafolds. 

Example 5.2. The Ladder Fractafold. Here T is the ladder graph consisting of two copies of 
Z, {a k } and {bk} with a k ~ b k and T consisting of three copies of Z, {x k+ i }, {w k }, {y k+ i} 



Figure 5.1. A part of the infinite Ladder Sierpinski fractafold. 



b-x b bi 
Figure 5.2. T graph for the Ladder Fractafold 
with w k joined to x k _i,x k+ i,y k _±, and y k+ ^, where x k+ \ is the edge [a k ,a k+ i], x y+ i is the 

X_3 X_X II X3 



2 2 2 2 




Figure 5.3. T graph for the Ladder Fractafold 
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edge [b k ,b k+ i] and w k is the edge [a k ,b k ]. 

It is easy to see that the spectrum of — Ap is [0,6], with the even functions <fg(a k ) = 
fe(b k ) = cos k6 or sink9, < 6 < n corresponding to A = 2 — 2 cos 6* in [0,4] and the odd 
functions ^e(a k ) = —ipo(b k ) = cos k9 or sin k0, < 6 < n corresponding to A = 4 — 2 cos 9 in 
[2,6]. 

These transfer to eigenfunctions of — A ro 

V?e(^fe+i ) = <Po(y k+ i) — cos(k + \)9 cos \B or sin(/c + \)9 cos \9 
(peiwk) = cos k6 or sinA;# and 
Mx k+ i) = -ipg(y k+ i) = cos(k + \)9 or sin(A; + \)9 

4>e(w k ) = 

with the same eigenvalues. It is also easy to see that there are no £ 2 (T ) eigenfunctions cor- 
responding to A = 6 (or for any A value whatsoever). Thus — Ar has absolutely continuous 
spectrum [0,6] with multiplicity 2 in [0,2] and [4,6] and multiplicity 4 in [2,4]. 

Example 5.3. The Honeycomb Fractafold. Here V is the hexagonal graph consisting of the 
triangular lattice L generated by (1,0) and (|, ^) and the displaced lattice £> + (|, ^). We 

denote by a(j, k) the points j(l, 0) + ^) of L and by b(j, k) the points a(j, k) + (|, ^) of 
the displaced lattice, with edges a(j, k) ~ b(j, k), a(j, k) ~ b(j — 1, fc) and a(j, fc) ~ 6(j, fc — 1). 
The eigenfunctions of — A r will have the form 




Figure 5.4. A part of the infinite periodic Sierpihski fractafold based on the 
hexagonal (honeycomb) lattice. 

k)) = e 2m ^ u+k ^ 
<Pu,v( b U> k )) = -fe 2ni( - ju+kv ^ 

where (u, v) G [0, 1] x [0, 1] and 7 depends on u, v. Let 1 + e 2nlu + e 2mv = re %9 in polar 
coordinates (so r and 9 are functions of u,v). Note that < r < 3. Then the eigenvalue 
equation requires r ) 2 = e 2t0 or 7 = ±e 2d with corresponding eigenvalues A = 3 =F r (so the 
choice ± yields the intervals [0,3] and [3,6] in spect(— Ar)). 

We can write the explicit spectral resolution as follows. For / e £ 2 (T) define 

/ a (v)^Je- 2 ^/(a(j,«)) 

3 k 
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6(-l,0) 



6(0,0) 



6(1,0) 



a(l,0) 



6(0,-1) 



6(1,-1) 



Figure 5.5. A part of the Hexagonal graph 



and 



fM = EE e_27ri °' u+M /( 6 o»)- 

j k 

We can invert these so that 

(f(a(j,k)y < - ! 

1/(60', *)) 



o Jo 



I 1 e ^i3u + kv)h faM + e -*» fb{UiV))dudv 
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+ 



e 2«o-«+fct;) -(f a ( u ,v) - e- i6 f b (u,v))dudv. 



Define A±(u, t>) by 

A± (u, f ) = 3 =F \/3 + 2cos27r-u + 2cos27ru + 2cos27r(w — v). 

For < A < 3 we define uq and i>0 by solving X + (u,v) = A, and similarly for 3 < A < 6 we 
solve \-(u,v) = A. We then define 



<51) WW,*))/"/. W"" 



27ri(j« 9 +fci; e )_ 



■(/o(t*e,Utf) ±e te f b (uo,v e )) 



d(u e ,v e ) 



d(x,e) 



de 



to obtain / = J Q 6 P\fd\ with —A r P x f = \P\f- This solves problem (a). 
To solve problem (b) we identify the space E 6 in £ 2 (r ). 




Figure 5.6. A part of the graph T for the Honeycomb Fractafold 

We may regard To as an infinite union of hexagons, each vertex belonging to exactly two 
hexagons. For any fixed hexagon H, define ipn to alternate values ±1 around the verticies 
of H, and to vanish elsewhere. It is easy to see that tpu in is E 6 . If {Hj} is an enumeration 
of all the hexagons in T then J2cj^ Hj (finite sum) is in E 6 . 
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Figure 5.7. Labels of hexagons and points 



Lemma 5.4. Suppose u 6 E e has compact support. Then u = ^CjipHj (finite sum). 
Proof. Suppose supp{u) C [^J Hj We will show that there exists jo £ A and c JO such that 

supp{u — Cj ipH n ) ^ [J ^j- The proof is then completed by induction. 

We choose j so that Hj lies in the top row and right-most down-right slanting diagonal 
of \Jj e AHj. In Figure [5771 above, j' = and u vanishes on Hi, H 2 , and H 3 . So u(x{) = 0, 
u{x 2 ) = 0, m(x 3 ) = 0. But u(x 3 ) + m(x 4 ) + u(y 3i ) = because E 6 = /cer(S' 2 ) and ^(7/34) = 
since 2/34 G f/3. So 14(2:4) = 0. A similar argument shows w(x 6 ) = 0. The only vertex left in 
H is £5. By subtracting off u(x^)ipH 5 we can make w vanish on if - 

We can systematically go across the top row in supp(u) from right to left and remove each 
hexagon, only changing u on the row below it. Eventually u will be supported on just one 




Figure 5.8. A row of hexagons 

row, and u(x) = unless x is one of the dotted points in Figure I5~gl 

Let H be the right most hexagon. The u \hi— implies u(xi) = and u{x Q ) = 0. 
Considering the triangle below the row we get u(x§) = 0. Considering the triangle above £4 
we get u(xi) = 0. So m \h = 0. □ 

Corollary 5.5. A function of compact support is in E§ if and only if u(x\) +u(x2) +u(x 3 ) = 
for every triangle {xi,x 2 ,x 3 } in Tq. 

Proof. The identity clearly holds for each ipjj, hence for all compactly supported functions 
in Eq. Conversely, every point x in r lies in exactly two triangles. Summing the identity 
for those two triangles yields the 6-eigenvalue equation at the point x. □ 
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The functions {ipH } do not form a tight frame, and it seems unlikely that they even form 
a frame (the lower frame bound is doubtful), so they do not seem well suited for describing 
Pq. We can, however, find an orthonormal basis of E 6 that consists of translates of a single 
function, but we pay the price that the function is not compactly supported. 

We change notation to index the hexagons in Figure l5\6l by the lattice [j, k] = j{i}+k{ ^ } . 

~2~ 

Note that hexagon Huu has six neighbors Hut ^n for 

If, k'] = [j, k] + {[1, 0], [-1, 0], [0, 1], [0, -1], [1, -1], [-1, 1]}. 
To describe a function 

(5.2) F = Y,f{{j,m H[m 

I? 

it suffices to give the discrete Fourier transform f(a, b) for (a, b) G [0, 1] x [0, 1] given by 

(5.3) f(a,b) = Y,f(lJ^])e- 2 ^ +bk \ 

for then 

(5.4) f([j, k})= I I e 2 ^ +b ^f(a, b)dadb. 



o J o 



In fact we will construct f(a,b) directly, and then substitute this in ( 15.41) and then in (15.21) 
to obtain our function in E e . 

The basic observation is that each point in Tq lies in exactly two neighboring hexagons, 
and the values of ipH for those two hexagons will be ±1. Thus 

<F,F>, 2(ro) = imw -mf, k'\)\ 2 

for / of the form (15 .2p . where the sum is over all neighboring pairs, and by polarization 

(5.5) < F,G >, 2(ro) = k}) - f([f, k'])(g([j, kj) - g([f, k 1 })) 
if F and G are of the form (15 .2p . Now we substitute (15 .4p in ( 15. 5 p to obtain 

(5.6) <F,G>, 2(ro) = 

[ I J2 e 27Ti{a]+bk) f(a, 6)[6-e 2 ™-e- 27ria -e 2mfe -e- 2 " 6 -e 27ri(a - 6) -e 2 " (b - a) ]^([j, k])dadb 
Jo Jo Z 2 

because of the form of the neighboring relation between [j, k] and [j', k']. But then we can 

evaluate the sum in (I5.6p using ( 15. 3p to obtain 

(5.7) 

1 /•! 



o Jo 



< F, G >£2( ro )= / / 2 (3— cos(27ra)— cos(27r6)— cos(27r(a — 6))) /(a, b)f(a, b)dadb. 



Lemma 5.6. The functions r Piq F = /([j, k] + [p, q])ipHy k] form an orthonormal basis of 

i? 

Eq for \p, q] e Z 2 if and only if 

(5.8) \f(a,b)\ = — 1 

y 2 (3 — cos(27ra) — cos(27r6) — cos(27r(a — b))) 
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Proof. We note that for T Ptq f(\j, k}) = f([j, k] + [p, q}) we have 

(5.9) (r p J)Xa,b)=e 2 ^ +b ^f(a,b) 
from (15.31) . so 

(5.10) < F,r Ptq F >£2 (ro) = 

f f e -2"(»P+6?)2(3_ C os(27ra)-cos(27r6)-cos(27r(a-6)))|/(a,6)| 2 dacZ6 
Jo Jo 

by (15. 9p and (15 .7L But the right side of (15. lQf) is 5(p, q) if and only if 

2 (3 - cos(2vra) - cos(2tt6) - cos(27r(a - b))) \f(a, b)\ 2 

is identically one, and this is equivalent to (15. 8 1) □ 

We are free to choose any phase in (I5.8p . but it is not clear what is to be gained, so we will 
simply choose f(a, b) to be positive. Note that the only singularity of / is near (0, 0), where 
it behaves like (a 2 + b 2 )~ 1 ^ 2 , so it is an integrable singularity, but not square integrable. Thus 
(15. 4 p is everywhere finite and decays like O ((j 2 + &; 2 ) -1 / 2 ). Although / is not in £ 2 (Z 2 ), we 
do have F e £ 2 (T ). 

Theorem 5.7. Let 

-1 pi e 2m(aj+bk) 



(5.11) f([j,k])= / / dadb. 

Jo Jo a/ 2 (3 — cos(27ra) — cos(27r6) — cos(27r(a — &))) 

Then < '^^T Ptq f([j,k])ipH [jk] ? is an orthonormal basis of Eg, and 

I Z2 J 

(5.12) 

w= E (E E T M hm)i>H aM (v)F{y)\ e ^M'^ H[f , k ^)- 

[p,q]ez 2 \yer [j,k]ez 2 J [j',k']ez 2 

Proof. This is an immediate consequence of Lemma 15.61 □ 

6. NON-FRACTAFOLD EXAMPLES 

Theorem 12.31 an be applied for examples that are not fractafolds. We assume that To = 
(Vo, E) is a finite or infinite graph which is a union of complete graphs of 3 vertices (it can be 
said that r is a 3- hyper-graph). In principle, we can allow T to have unbounded degrees, 
as well as loops and multiple edges, but in this section we will keep everything simple and 
assume that r is a regular graph. As before, each of these complete 3-graphs we call a cell, 
or a 0-cell, of r . We denote the discrete Laplacian on T by A ro . We define a finitely 
ramified Sierpinski fractal field $ by replacing each cell of T by a copy of SG. These copies 
we call cells, or 0-cells, of Naturally, the corners of the copies of the Sierpinski gasket 
SG are identified with the vertices of IV See [12] for fractal fields, not necessarily finitely 
ramified. Since the pairwise intersections of the cells of $ are finite, we can consider the 
natural measure on 5", which we also denote \l. Furthermore, since A^^ is a local operator, 
we can define a local Laplacian A on ^, in the same way as explained in [37] (this means 
that the sum of normal derivatives is zero at every junction points). One can see that most 
of our results can be easily generalized for the finitely ramified Sierpinski fractal fields. For 
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Figure 6.1. A part of the periodic triangular lattice finitely ramified 
Sierpihski fractal field. This fractal field is not a fractafold. 




Figure 6.2. A part of the infinite triangular lattice, the T graph for the 
fractal field in Figure 16.11 

instance, Theorem 12 .31 is essentially still valid. One change to be made is that on the graph T 
we have to consider the probabilistic Laplacian (which is explained in [2U| I3"3"]). and multiply 
it by 4 to align with the normalization of the Laplacian on the Sierpihski gasket. 

In the example shown in Figure IST21 the spectrum on To is [0, 8] for the adjacency matrix 
Laplacian, and the spectrum is [0,4/3] for the probabilistic Laplacian. Thus So = [0, -y]. 
In this particular case the spectrum is absolutely continuous by the classical theory (see 
[21] |22| l23~t |24"| |25| and references therein] for a sample of relevant recent results on peri- 
odic Laplacians). Combining the methods described in this paper, we obtain the following 
proposition (see also Figure [673]) . 




Figure 6.3. Computation of the spectrum on the triangular lattice finitely 
ramified Sierpihski fractal field. 



Proposition 6.1. The Laplacian on the periodic triangular lattice finitely ramified Sierpinski 
fractal field consists of absolutely continuous spectrum and pure point spectrum. The abso- 
lutely continuous spectrum is $H _1 [0, y]. The pure point spectrum consists of two infinite 
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series of eigenvalues of infinite multiplicity. The series 5fH _1 {3} C 9^ _1 {6} consists of iso- 
lated eigenvalues, and the series 5£R: 1 {5} = 9i -1 {0}\{0} is at the gap edges of the a.c. 
spectrum. The eigenfunction with compact support are complete in the p.p. spectrum. The 
spectral resolution is given by (I2.14p . 

It is straightforward to generalize such a result for other finitely ramified Sierpihski fractal 
fields (see, in particular, Remark [5 . 1 [) . 
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